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We have found that the conduction in Si-MOS structures has a substantial imaginary component 
in the metallic phase for the density range 6 x Uc > n > n^, where Uc is the critical density of 
the metal-insulator transition. For high mobility samples, the corresponding delay (or advance) 
time, r ~ (0.1 — 10) ms, increases exponentially as density and temperature decrease. In very low 
mobility samples, at temperature of 0.3K, the time-lag in establishing the equilibrium resistance 
reaches hundreds of seconds. The delay (advance) times are approximately 10^ — 10* times larger 
than the overall _RC-time of the gated structure. These results give evidence for a non-Boltzmann 
character of the transport in the low-density metallic phase. We relate the time-dependent effects 
to tunneling of carries between the 2D bulk and localized states. 



The metal-insulator transition observed in different 
two dimensional carrier systems jl) is currently in the 
focus of interest. A number of models put forward for 
its explanation span from a non-Fermi-liquid state 
1^ to interface traps physics In the latter mod- 

els, the strong exponential drop in the resistivity at low 
temperature is a transient temperature effect only. On 
the experimental side, from measurements at high car- 
rier density the exponential drop was proven not to 
be related to the ground state conduction, at least for 
high carrier densities n > (10 — 15) x ric- 

The involvement of interface traps into the transport 
may be revealed by studying charging effects, time lag, 
noise character etc. We present here data evidencing 
that the time-dependent effects are essential in the metal- 
lic phase, even for densities 6 times the critical one, 
Uc. We studied in detail four samples with different 
peak mobility, Si-11 (/LtP'="'= = 39, OOOcmVVs), Si-22 
(^peafe ^ 33,000cmVVs), Si-4/32 (/i = 8400cmVVs), 
and Si-52 (/i — 1300cm^/Vs). For the first three sam- 
ples, an Al gate film was deposited onto the Si02 layer 
followed by a post-metallization anneal. The density of 
trapped carriers (estimated from the threshold voltage 
0), was (2 — 3) X lO^^cm"^ for the first two samples, 
and 21 x 10^" cm^^ for the third one. In the lowest mo- 
bility sample, we increased intentionally the amount of 
disorder by thermal evaporating an Inconel gate without 
a subsequent anneal. Beyond the substantial decrease in 
the mobility, the amount of trapped carriers increased up 
to 60 X 10^" cm^^. All sarnples were of the same Hall- 
bar geometry, 5 x 0.8 mm^, 0| with gate oxide thickness 
of rfsi02 — 200 ± 20 nm, aspect ratio of w/l — 0.32, and 
corresponding capacitance between the gate and 2D layer 
C « 690 pF. The potential and current contacts to the 2D 
channel were lithographically defined and made by ther- 
mal diffusion of phosphorus. The overall device i?C-time, 
including contact resistance was of the order of (1—10) ^s, 
and was expected to contribute a negligibly small imag- 
inary component to the sample ac-conductance. 



Four-terminal ac-transport measurements were carried 
out in the frequency range 0.3 to 30 Hz with a quadrature 
lock-in amplifier. In order to eliminate the influence of 
the resistance of potential probes, we used a battery oper- 
ated electrometric preamplifier with an input current less 
than 1 pA. The amplifier phase-frequency characteristic 
was verified not to contribute to the studied effects. In 
all samples, we found the time-dependent effects to per- 
sist in the metallic phase, far above the critical density. 
In high mobility samples, the characteristic times were 
in the ms-range and were measured from a phase shift (p 
between the voltage drop, Vx, and the source-drain cur- 
rent, Ix, as well as from the frequency dependence f{F). 
In low mobility samples, the characteristic times were of 
order of 1-100 s and were measured directly, by applying 
a small voltage step on the top of the constant gate volt- 
age and measuring the transient voltage Vx at a constant 
current ■ In the following, we define the "resistivity" , 
p as the in-phase component, Re{{Vx/Ix) x w/l). 

Measurements with high mobihty samples. Fig- 
ure 1 shows the phase shift between the ac voltage Vx 
and the current Ix, measured in a high mobility sample 
at a frequency of 3.8 Hz, as a function of carrier density. 
The phase shift emerges as the density decreases below 
6 X 10^^ cm~^. This is about 6 times larger than the criti- 
cal density, ric = 0.95 x 10^^ cm~^ for the metal-insulator 
transition in this sample. 

The lower inset of Fig. la shows that, as density de- 
creases, the phase shift first is negative (which corre- 
sponds to the voltage delay), then becomes positive (volt- 
age advancing), and, finally, becomes again negative close 
to the critical density. The phase shift increases lin- 
early with ac-current frequency F and may thus be in- 
terpreted as a time delay (or advance, correspondingly), 
T = (p/{2TrF). The upper inset in Fig. la shows delay 
time {— \t\) calculated from the slope of the frequency 
characteristics, at a fixed temperature of 290 mK, and 
over the range of high densities where (f) < 0. As temper- 
ature decreases, the phase shift displays more and more 
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pronounced oscillations as a function of density. The 
phase shift is reproducible during the same cool-down, 
however, in different cool-downs the oscillatory details 
varied on the density scale. 
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FIG. 1. (a) Phase shift between ac- voltage Vx and 
current Ix measured with sample Si-22 at frequency 
3.8 Hz for 6 temperatures (indicated on the figure). 
The range of n = (4 — 7) x 10^^ cm~^ is blown up 
on the lower inset. Upper inset displays delay time 
vs density, measured from the frequency dependence 
of the phase shift at T = 0.29 K. (b) "Activation en- 
ergy", To, for low and high temperature ranges. Ver- 
tical dashed lines n = rii separate the regions of pos- 
itive and negative ip. 

Figures 2 a and 2b show the phase shift and the resis- 
tivity for sample Si-22 as a function of temperature for 
eight fixed densities. It is remarkable that the strong ex- 
ponential drop in resistivity develops in the same ranges 
of densities and temperatures as the phase shift does, al- 
though we can not simply relate the two effects to each 
other. 

As follows from Figs, la and 2a, both, r and ip decay 
about exponentially with density and temperature. 



T oc /i(n,r) exp(To(n)/T). 



(1) 



The prefactor /i oscillates as a function of density (and 
of temperature) changing sign at "node" values, n^. The 
definition of the slope, Tq, is illustrated by the dashed 
tangent lines in the lower inset of Fig. 2a. 

TQ{n) is not constant over entire temperature range: it 
is large for high temperatures, and decreases for lower 
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FIG. 2. (a) Temperature dependence of the phase 
shift, <p, measured at 3.8 Hz for eight different den- 
sities (indicated on the figure in units of lO'^^ cm~^). 
Lower and upper insets are the "Arrhenius plots" for 
<p vs 1/T: lower inset is for the negative ip over the 
temperature range "B" ; upper inset is for the positive 
ip over the range "C". (b) Temperature dependence 
of the resistivity for the same eight densities as in (a) . 



T's. In the vicinity of the nodes ni and for low tem- 
peratures, the slope tends to vanish which is simply af- 
fected by the oscillatory behavior of fi{T). Therefore, 
the narrow density ranges around the nodes were ignored 
in the calculations of Tq. The resulting density depen- 
dence To(n) is shown in Fig. lb, evaluated separately for 
high (T > 2.5 K) and low (T < 2.5 K) temperatures. De- 
spite Eq. (1) describes the data very roughly, an impor- 
tant conclusion can be drawn immediately: Tq docs not 
decrease to for the nodes n — n.i and develops smoothly 
from the ranges of > to those oi tp < 0. This means 
that the nodes are related to the prefactor fi{n) rather 
than to the exponential factor. For high densities, Tq 
decays steeper than oc (n — uq)^^ , as shown in Fig. 1 b, 
thus causing an exponential decay of r for high densities. 

Low mobility sample. In the low mobility sample 
Si-52, the time-dependent effects are much stronger and 
manifest themselves in a transient voltage between poten- 
tial probes when the gate voltage Vg changes by a small 
step AVg <C Vg. Typical transient curves i - 5 of Ap{t) = 
p{t) — p(0) normalized by Apo = p(oo) — p(0) are shown 
in Fig. 3a, for 5 different densities. The curves were fit- 
ted with exponential functions po exp(— t/r^) (shown by 
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dashed curves), from which the time lag Td was obtained. 
At some, rather arbitrary densities, the transient curves 
were non-monotonic with oscillations (curve 5), or jumps 
(curves 7, 8, 9). 
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FIG. 3. (a) Ap(t), normalized by p(oo) - p(0). 
Curves 1-5 are shifted vertically for clarity. The 
curves 1-5 correspond to the carrier densities 3.3, 
5.5, 7.7, 12.1, 23.1 (in units of lO^cni"^). Curves 
6-9 {n = 34, 26, 23) in the inset show examples of 
non-monotonic behavior. (b) Time lag (right 
y-axis) and p (left j/-axis) vs carrier density. 

The time lag and resistivity p are plotted in Fig. 3b 
vs carrier density. Their ratio, r^/p, is again 10* times 
larger than the sample capacitance pointing to its irrele- 
vance. 

Discussion. The characteristic times, even for high 
mobility samples, are of the order of 10~^ — 10~^s and, 
therefore, can not be associated with any RC time con- 
stants of the sample. Having the typical sample parame- 
ters R - (10^ - 10^) Ohm, and C = 0.7 x lO^^ F, one can 
hardly find in the sample either a capacitance ~ 10~* F, 
or a resistance ~ 10* Ohm (in the metallic range of den- 
sities). Only for high densities n > 20 x lO-'^^cm^^, the 
delay time becomes comparable to RC = 10~^ s. The 
huge time r can not be attributed to contact phenom- 
ena, because for the more disordered sample Si-52, at 
much higher densities (and for lower resistance of the 
contacts, correspondingly), r is larger by a factor of 10^. 
The irrelevance of the contacts to the time-dependent ef- 
fects is also confirmed by the linearity of the / — ^-curves 
measured between different contacts with currents in the 
range from 10^^ down to 10^^^ A. 

Model. Interface defect charges originating from the 
lack of stoichiometry are intrinsic to Si/Si02 system; 
their typical density is lO^^cm"^ for a state-of-art ther- 
mally grown dioxide . Tunneling of electrons from Si to 
the interface charged states is known to cause a time lag 
in Si-MOS capacitors at room temperature These 
charged states are partly neutralized during slow cool- 
ing of Si-MOSFETs with a positive gate voltage applied. 



Further, at liquid helium temperatures, electron tunnel- 
ing rate to these interface traps in Si02 under a barrier 
of 3.2 cV is negligible. The uniform part of the potential 
produced by interface-state charge is out of importance, 
however, spatial fluctuations of the built-in charge pro- 
duce shallow fluctuations of the potential acting on 2D 
electrons in Si (at z > 0), and cause corresponding local- 
ization of electrons. We assume that the observed times 
are due to tunneling processes, on the Si-side entirely, be- 
tween the electrons in 2D "bulk" and the potential traps 
in the localized areas produced by the fluctuations of the 
interface charge For the discussed low-temperature 
case, the active traps are created by the attractive (pos- 
itive) charges that fall inside a large-scale repulsive fluc- 
tuation. The attractive charge placed at the interface 
(at z < 0) localizes an electron nearby (on the Si-side, 
z > 0), with a binding energy |^ 

Eb = -m*e^/8K^n^. (2) 

Here to* ~ 0.21me is the electron effective mass, k — 7.7 
the average dielectric permittivity, and thus Eb — 0.02 eV 
for the Si/Si02 interface. The repulsive charges which are 
very closely located to the attractive charge, decrease the 
binding electron energy. As a result, the binding energy 
distribution broadens and extends over a wide energy 
range, from about Eb down to 0. The effective binding 
energy e^^'^ is therefore substantially lower than Eb- 

The localized state is located inside a large-scale repul- 
sive fluctuation and is surrounded by a broad potential 
barrier of the height ejj'^^. The barrier itself is surrounded 
by the electrons in the metallic regions of the 2D bulk. 
The barrier is responsible for the large electron capture 
and emission times. At nonzero T, only the traps located 
close to the Fermi level, within Ep ± kT, are recharging 
when the local potential in the 2D bulk varies. 

For low temperatures, the electron emission time |^ 
equals to 

Tem « (?i/£6)exp(-a;/A), (3) 

where A = •y/s^/8m*e^^ is the typical tunneling length 
and X is the distance from the trap to the nearest con- 
ductive region. The capture time Tc is related to T^m 
by the obvious relationship: Tc/Tem = (1 ~ f)/fj where 
/ is the level occupancy (Fermi distribution function). 
For the traps whose energy is close to Ep, f 1/2 and 
these two times are about equal. Using, for an estimate, 
^b'^'^ = O.OleV we obtain A — 20A and the tunneling 
distance x — 460A corresponding to the tunneling time 
10"'^ s. Thus, a radius of the repulsive barrier r is to 
be of the order of 500A, to account for the recharging 
time of 1 ms. For more disordered samples, the amount 
of the charge trapped at the interface and the amplitude 
of potential fluctuations are even larger. Therefore, the 
radius (in the x — y plane) of the potential fluctuations 
is also larger. To account for r = 100s, we estimate x 
has to be equal to TOOA. The length scale, (500 - 700)A, 



3 



of the potential fluctuation is consistent with numerous 
data obtained for similar samples, as well as with direct 
tunneling microscopy of the Si/Si02 interface [ p^ . 

As electron density and Ep increase, the potential bar- 
riers get thinner due to screening by free electrons of 
the 2D bulk, and the radius of the total repulsive fluc- 
tuation decreases, causing r to decrease monotonically. 
This is in a competition with a density dependence of 
the effective binding energy el^'^ ; as a result, the overall 
density dependence of t may be non-monotonic. Finally, 
above a certain density, all localized states sink below 
the Fermi energy, the barriers are screened entirely, and 
the Drude-Boltzmann regime sets in replacing tunneling. 
The onset occurs at Fermi energy equal roughly to the 
amplitude of bare fluctuations, and is thus inversely pro- 
portional to the sample peak mobility. 

The temperature comes into the model via (i) the tun- 
neling distance to the nearest trap level found within the 
energy interval Ep ± kT, (ii) Fermi distribution function 
/ and (iii) activation processes on the border and inside 
the localized area. These mechanisms lead to tempera- 
ture dependences r oc exp(— To,i/T), and the resulting 
temperature dependence may have different Tg^i in dif- 
ferent ranges of temperature. 

The electron tunneling time is much larger than the 
transport scattering time 3ps), the electron-electron 
interaction time, h/Eee ~ 0.1 ps, and the electron diffu- 
sion time (« 20 ps). Therefore, all electrons in 2D layer 
do participate in tunneling during a time r > 10~^ s. The 
capacitance and Hall voltage are measured at fre- 
quencies lower than the tunneling rate, 1/rem, hence, all 
the electrons of the 2D bulk participate in re-charging or 
in Hall transport. 

In summary, our data show that the charge transport 
in the "metallic conduction" regime is accompanied by 
(or includes) a non-diffusive component. It manifests it- 
self in the time-dependent effects in metallic phase over 
density range from Uc to about 6 times Uc. The de- 
lay/advance time between the ac voltage and the current 
in high mobility samples is of the ms-range and grows 
exponentially as density and temperature decrease. For 
more disordered samples the time lag between the gate 
voltage pulse and the response reaches 1-100 seconds. We 
associate these times with the "in-plane" tunneling of 
carriers between the 2D bulk and the potential traps. 
The suggested model presumes that the unit "slow" trap 
consists of a potential well surrounded by potential bar- 
rier, and seems to be applicable to various material sys- 
tems since similar long-range localized areas were found 
in GaAs/AlGaAs as well ||l^. Particularly, this may be 
relevant to the system with a set of artificial quantum 
dots (traps) ||l^]. The model may qualitatively describe 
(i) large delay time r, (ii) the growth of r as tempera- 
ture and density decrease, and (iii) the non-monotonic 
density and temperature dependence of t. Whereas the 
tunneling time in the above model is different from that 
in Ref. 0], the physics of the temperature dependence 



of the resistivity caused by "fast traps" (located nearby 
the border between the localized and free carriers) may 
be similar. The complexity of the presented data, how- 
ever, requires a thorough theoretical consideration, which 
should take into account interaction of the free carriers in 
2D bulk with the localized ones, recharging and screen- 
ing processes within the localized areas, and symmetry 
properties of the surface localized states. 
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